Solitons in Poly(oxyethylene) 
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We introduce a soliton model to describe conformational structure of poly(oxyethylene) in aqueous 
solution. The model is based on former investigations for twistons in polyethylene, and the soliton 
solution corresponds to stable excitation that locally transforms the gauche conformation of the 
O-CH2-CH2-O dihedral angle from the positive (negative) configuration to its negative (positive) 
partner. The conformational deformation in the polymer main chain is mediated by the unstable 
cis conformation which changes the physical profile of the chain, contributing to trigger important 
properties of the polymer. For instance, it may activate the process for the capture of metallic ions 
in aqueous solution and it may also enlarge solubility of the polymer in water. The soliton solution is 
used to guide a quantum chemistry investigation that examine, and characterize, a chain containing 
40 monomeric units, tracing quantitative correlations with the soliton model. 

PACS numbers: 61.41. +e, 03.50.-z, 05.45.-a 



Solitons appear under the direct influence of nonlin- 
earity, and may play important role in several branches 
of nonlinear science. In particular, they may spring as 
finite energy excitations in nonlinear polymeric chains. 
A classic example is polyacethylene (PA): in this chain, 
single and double bonds alternate, generating Peierls in- 
stability that gives rise to solitons, as accounted for in 
the Su-SchriefFer-Heeger model [Q; see Ref. |^ for fur- 
ther details. As one knows, the PA chain admits an in- 
teresting field theoretical formulation in terms of bosonic 
and fermionic fields ^, ^. Another example of a nonlin- 
ear polymeric chain that supports solitonic excitations is 
polyethylene (PE): in this case the chain only contains 
single sigma bonds, which activate degrees of freedom of 
the torsional type, giving rise to twiston excitations in 
the polymeric chain . The presence of twistons in the 
PE chain was further explored in Refs. |^], and more 
recently in Ref. |^ one has shown that these twistons can 
be seen as topological solitons that appear in a specific 
model of two real scalar fields. 

The work of Ref. [|| follows the investigations on topo- 
logical solitons iniciated in |^ , and further explored more 
recently in |l^, |l2j . The main idea is that the pres- 
ence of topological solitons in a polymeric chain may di- 
rectly contribute to the conformational characteristics of 
the polymer, a fact important to the understanding of 
physical, chemical and biological properties of the sub- 
stance. A singular example is the poly(oxyethylene) 
(POE), which has long been the subject of great atten- 
tion due to its various biotechnical and biomedical appli- 
cations In this letter we develop the idea that the 
POE chain is able to support topological solitons, and 
that these solitons contribute significantly to the under- 
tanding of important features of this polymer. 

We start with the PE chain, which is a very long chain 
made of CH2 units. It may be seen as the limit for large 
n of the chain family [(CH2)n-0-]m. In this family the 



second member (n = 2) is the POE, which is formed by 
repetitions of the chain unit CII2-CH2-O-. The POE is a 
very important polymer in this family, because the pecu- 
liar array in the chain provides to this polymer some very 
impressive features, among them the remarkable ability 
to entrap metallic ions in aqueous solution and its solu- 
bility in water to almost any extent. Besides, one knows 
that ethylene oxide polymers are amongst the most com- 
monly used substances in pharmaceutical and other in- 
dustrial formulations ||, 0, |l|, |l|, 0, |20| . 

Molecular dynamics simulations [|l^ |l|, ^ on 
the conformational structure of the POE and 1,2- 
dimethoxyethane (DME = CH3-O-CH2-CH2-O-CH3, 
which is correlated with the POE's dimer) have shown 
that the polymeric chain may present a very rich fam- 
ily of conformational possibilities, with preference for tgt 
conformers in condensed media - we are using standard 
notation, where 't', 'c' and '5' represent trans (tt), cis 
(0) and gauche (tt/S) arrangements, respectively. The 
tgt conformer was also observed in POE chains used 
as hidrophilic moiety of nonionic surfactants, in aque- 
ous solutions, through Raman spectroscopy |23|. The 
stabilization of the tgt conformer in the liquid or solid 
state may be understood in terms of polar intermolec- 
ular interactions, because this conformer has an appre- 
ciable dipole moment while the ttt conformer (all trans 
conformations) has zero dipole moment [T^ . This pref- 
erence for tgt conformation is also predicted for DME by 
experimental NMR |^ and infrared spectroscopic 
studies. The water affinity of POE is directly related to 
a structural similarity between the POE chain and liquid 
water, which is not found for other polyethers p7| : 
when substantially hydrated, the POE chain assumes lo- 
cal conformations mostly with the tgt arrangement; in 
this case there is almost perfect match between the dis- 
tance of the oxygen atoms in the structure of liquid water 
and the distance between the nearest neighbor oxygen 
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atoms in the polymer chain. Thus, the POE chain may 
be easily incorporated into the structure of liquid water, 
with the ether oxygen atoms replacing some of the water 
molecules, with the hydrophobic ethylene groups being 
encaged in voids in between water molecules. There are 
other conformational possibilities in the POE chain; for 
instance, the tgg conformer is considered the second most 
populated in DME, also based on spectroscopic results 
both for DME and POE, probably because of the pres- 
ence of an important intramolecular 1,5-CH- • -O inter- 
action that stabilizes this particular conformation. The 
nomenclature ttt, tgt, and so forth, that we are using is 
valid to represent the three most important dihedral an- 
gles in the polymer main chain, namely O-CH2-CH2-O, 
CH2-CH2-O-CH2, and CH2-O-CH2-CH2, as we illustrate 
in Fig. 1. 
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FIG. 2: The cis conformation that intermediates the most 
stable gauche conformation for the O-CH2-CH2-O dihedral 
angle. The central circles are drawn to represent the two 
neighboring carbon atoms between subsequent oxygen atoms. 




FIG. 1: The three most important dihedral angles in the POE 
main chain. 

We can visualize topological solitons in the POE chain 
thinking of the tgt conformation, the most stable confor- 
mational configuration in the POE chain in liquid water. 
The symmetry of the chain admits two equivalent possi- 
bilities, where the helix can be clockwise, generated by 
a sequence of tg^t arrangements, or counterclockwise, 
generated by a sequence of tg~t orientations. These two 
gauche conformations (helixes) are degenerate in energy, 
and they are typically very populated, according to the 
previous discussion 1 24[ 2^ . They appear to be the most 
stable of all the conformational possibilities the POE 
chain may comprise. Thus, although the POE chain may 
adopt other distinct conformations, our model considers 
only the two most important set of degrees of freedom, 
represented by the O-CH2-CH2-O and CH2-CH2-O-CH2 
dihedral angles, and the most stable conformation, the 
gauche and trans arrangements, respectively. For both 
angles, the cis conformation is unstable, and is seen as 
a barrier for gauche and trans orientations to tunnel to 
its symmetric partner. We illustrate this in Fig. 2 for 
the gauche conformation of the O-CH2-CH2-O dihedral 
angle. 

This reasoning allows introducing a model of two real 
scalar fields to describe the nonlinear chain. Under this 
assumption we could establish, a priori, the two most rel- 
evant sets of degrees of freedom to describe the transition 



tg~^t ^ tg~t that springs as the topological soliton. 

The connection between these detailed information at 
molecular level with the mathematical representation of 
the polymer is implemented by the field theoretical model 
that we now develop. The model is described by two real 
scalar fields, (f) and x, which live in the bidimensional 
space-time {x,t). The Lagrangian density has the form 

/: = ^5„(/.9"0-f i9„x5"x-mx) (1) 



where 



(2) 



Here da(j) — d(j)/dx", ~ dW/dcfi, etc. We are using 
standard notation, with = {x'^ = t,x^ — x), and 
Xa — {xq = t,xi — —x), and ?i = c = 1. The form of 
the potential is special, because it allows obtaining static 
solutions of the equations of motion 



= w^w^^ + w^w^^ 
= w^w^^ + w^w^^ 



through the first order differential equations 



dx 
dx 
dx 



(3) 
(4) 

(5) 
(6) 



Solutions of the first order equations are stable and 
are minimum energy solutions - see ]T]| and references 
therein for more details. 

The model we consider is defined by the superpotential 
W = ]¥{(/), x), which has the specific form 



(7) 
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We consider A, /x, v real, and a real and positive. This 
model is new, and reproduces the model first considered 
in Ref. lotl in the limit ^ 0. We calculate Wa,a, and 



to see that 



(8) 



We notice that the superpotential is harmonic for ^+A = 
0, thus the equations of motion can be reduced to a family 
of first order equations for A + ^ = |^ . 

The potential presents absolute minima at the points 
where Xi/o = va/2iJL(j)i are given in terms of 4>i, 
which represent the real solutions of 



- - h — 



(9) 



There are two solutions for i^^/Xfi < 0, four solutions 
for < ly^/fiX < 1, two solutions for v^/\^ — 1, and 
no solution for j [i\ > 1. We consider the case < 
ly^/fiX < 1. We write /i = A?', and — As, to get to 
< s^/r < 1. This implies that r > 0, and we take 
-~^/r < s < V^. 

The equations of motion for static fields are solved by 
solutions to the first order differential equations, after 
changing (f) ^ a(j), x ^ o^Xi x Xax 



d(j)_ 
dx 
dx 
dx 



1 - 0^ - rx^ 



— = s — 2r(j)x 



(10) 
(11) 



We notice that this system does not allow one-field solu- 
tions for s 7^ 0. This peculiar behavior of the model is 
important, and is in good agreement with the quantum 
chemistry investigation that we have implemented, as we 
comment on below. We can find true two-field solutions 
to the above first order equations. They are described by 
the field configurations 



(j)s(x) 

Xs{x) 



a\l — tanh(Aav2rx) 
2r 



(12) 



±a 



2r 



2 tanh(Aa\/2rx) (13) 



which are valid for s — ±-\/2 — l/r, with r > 1/2. We no- 
tice that the limit r 1/2 transforms the above solution 
back to the one-field solution of the simpler model, where 
s = 0. We see that the amplitude of (j) is — a/V2r, 
and that of x reads = |s | A^. Thus, for 1/2 < r < 1, 
we get < < 1, and the amplitude of x is always 
smaller then that of </>. For r > 1 we get > 1, and 
in this case the amplitude of x is greater then that of 4>. 
A complete investigation of the two-field model will be 
published elsewhere. 

In the POE polymer, the topological soliton represents 
a conformational distortion in the tgt chain: it describes 



the change from the tg~^t conformation to the tg~t con- 
formation. This modification is done with the chain pass- 
ing through the unstable cis conformations, in a way such 
that the topological distortions locally increase the dipole 
moment in the chain, which are evenly distributed along 
the entire chain in both the tg~^t and tg^t conforma- 
tions. The presence of a topological soliton breaks the 
dipole moment uniformity in the chain, increasing the 
polar nature of the chain in the soliton core, and this 
is the mechanism that triggers the process for captur- 
ing metallic ions in aqueous solution, or that favors the 
presence of hidrogen bonds in water solutions. 

To make our investigations more quantitative we now 
examine the connections between the detailed informa- 
tions obtained in the model described by two real scalar 
fields and the POE chain at molecular level. This can be 
obtained through quantum chemistry calculations, with 
the mathematical representation of the polymer as given 
by the field theoretical model. We investigate the chain 
H(-CH2-CH2-0-)4oH, a POE chain with 40 monomeric 
units. We suppose that the topological soliton repre- 
sents the transition state between the positive and neg- 
ative helix conformations for the polymer main chain. 
The transition state between the two most stable con- 
formations of the POE is represented by a polymer that 
mixes the gauche and the energetically unfavorable cis 
conformations at the same molecule, in a way such that 
it makes possible the application of quantum chemistry 
machinery to do the investigation, since the soliton ex- 
citation is a continuum phenomenon that needs an ex- 
tended molecular representation. In Fig. 3 one depicts 
the polymer molecules considered as reagent, named pos- 
itive helix (lower figure), and as transition state, named 
soliton (upper figure). We established that the dihedral 
angle C-O-C-C was maintained at configuration trans (tt) 
throughout the entire polymer molecule, and that the 
other two dihedral angles in the main chain - namely 
O-CH2-CH2-O and CH2-CH2-O-CH2 - are choosen ac- 
cording to the field theoretical model, in which we are 
identifying the cf) and x fields with the two dihedral an- 
gles CH2-CH2-O-CH2 and O-CH2-CH2-O, respectively. 
[In fact, we have done several other investigations: for 
instance, if the CII2-CH2-O-CH2 angle is also frozen, the 
results are very unusual, showing that the polymer main 
chain visit very energetic configurations, returning along 
the same chain axis, but reversing its direction; also, if all 
the three dihedral angles are free, the polymer main chain 
bends considerably, breaking the unidirecional character 
of the polymer] . 

We use results of the two field model to see that the 
modification of the two dihedral angles that are corre- 
lated with important degrees of freedom of the main 
chain follow the hyperbolic tangent behavior. We identify 
the two dihedral configurations O-CH2-CH2-O and CII2- 
CH2-O-CH2 with the X and (p fields, respectively, and 
so we use the quantitative pattern of the soliton solution 
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FIG. 3: The two conformational structures for the POE main chain. The lower chain is shown in the positive helix conformation. 
The upper chain shows the solitonic structure, which drives the most important dihedral angles in the POE main chain, changing 
the gauche conformation from positive (negative) to negative (positive) helicity. 



( p^ and (|T^) to set the values for dihedral angles through 
the 40 monomeric unities of the POE molecule. The posi- 
tive helicity is achieved by setting the O-CH2-CH2-O and 
CH2-CH2-O-CH2 angles with the values -|-7r/3 and — tt, 
respectively, in the entire polymer molecule - see Fig. 3, 
lower chain. For the soliton, however, we consider the 
first 8 monomeric units in both extremities of the poly- 
mer as positive and negative helicities, respectively, while 
for the 24 central units the values of the O-CH2-CH2-O 
and CH2-CH2-O-CH2 dihedral angles vary according to 
the hyperbolic tangent profile. The O-CH2-CH2-O angle 
was modified from -I-tt/S to — 7r/3, while the CH2-CH2- 
O-CH2 angles was modified from — tt to +7r, through the 
curvature of each one of the two hyperbolic tangent so- 
lutions that describe the soliton. This identification re- 
strains the parameters r, s, and a that appear in the field 
theoretical model to the values r = 9/17, s = —1/3, and 
a = TT yJl'St/n and the soliton solution can be rewritten 
as 



TT tanh ( a/18/17 xj (14) 
X{x) = -(7r/3) tanh(^Vl8/17 (15) 

where we have changed \ax x, to use dimensionless 
coordinate. 

The molecules shown in Fig. 3 represent the final re- 
sults obtained with semi-empirical molecular orbital cal- 
culations that were made by using AMI |^ method im- 
plemented in the MOPAC 93 version 6.0 package 
During the geometry optimization procedures used in or- 
der to find the stationary point at the potential energy 
surface (PES) in each case (positive helicity and soliton), 
all the dihedral angles were maintained frozen, while the 
normal angles and bond lengths were allowed to relax. In 
Fig. 3 we notice that the upper chain (which contains the 
soliton) slightly deviates from its unidirensional charac- 
ter, but the linear configuration is still in good approx- 
imation with the field model. The main results for the 
electronic properties are given in the Table below, which 
shows the heat of formation and dipole moment (fi) for 
both the soliton (upper line) and the helix (lower line) 
molecules. 



Heat of Fomation (kJ) 


fi (Debye) 


-7621.47 


11.53 


-7719.56 


1.42 



The results reveal that the soliton is 98.09 kJ more en- 
ergetic than the normal gauche conformation, in agree- 
ment with the fact that the soliton induces a localized 
transition between positive and negative helicities of the 
gauche conformation. This energy difference is explained 
by the instability of the cis or near-cis conformations in 
both O-CH2-CH2-O and CH2-CH2-O-CH2 dihedral an- 
gles, found in the soliton core. We believe that it can 
be minimized in condensed media, particularly in polar 
solvents like water, for example, because of the greater 
polarity of the soliton molecule (/i = 11.53 D) in com- 
parison to the positive or negative helicity (/i = 1.42 D). 

The soliton model that we have used seems to be lim- 
ited, since it does not account for intermolecular inter- 
actions between neigbour POE chains. However, this is 
a good approximation in aqueous solutions, even at high 
concentrations, since interactions between POE chains 
via hydrogen bonds are scarce In fact, there are 

more interactions between oxyethylene units in the same 
chain ||2^, mediated by hydrogen bonds via neighbour 
water molecules, but this we have empirically taken into 
account by considering the positive and negative helix 
conformations for the POE main chain as the most pop- 
ulated conformations in aqueous solutions. 
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